# p02_dem_resamp.sh
# - 描述: 读取 DEM 并进行重采样
# - 作者: iam002
# - 日期: 2025年2月24日
# - 版本: V1.0

# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
# 读取全局变量
# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
if [ ! -f "$1" ]; then
    source p00_global_params.sh
else
    source $1
fi

# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
# 日志路径
# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
cur_script_name=$(basename "$0" .sh)
out_log_path=$in_output_folder/$cur_script_name".log"

# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
# 设置geoid
# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
if [ "$in_use_geoid" == "egm2008" ]; then
    geoid_dem_path="$DIFF_HOME/scripts/egm2008-5.dem"
    geoid_dem_par_path="$DIFF_HOME/scripts/egm2008-5.dem_par"
elif [ "$in_use_geoid" == "egm96" ]; then
    geoid_dem_path="$DIFF_HOME/scripts/egm96.dem"
    geoid_dem_par_path="$DIFF_HOME/scripts/egm96.dem_par"
else
    geoid_dem_path="-"
    geoid_dem_par_path="-"
fi

# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
# 读取 DEM
# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

tmp_dem_path=$in_output_folder/$in_produce_name".dem"
tmp_dem_par_path=$in_output_folder/$in_produce_name".dem_par"


echo "[input] DEM tiff path: " $in_dem_path | tee $out_log_path
echo "[input] DEM resample factor x: " $in_dem_resamp_factor_x | tee -a $out_log_path
echo "[input] DEM resample factor y: " $in_dem_resamp_factor_y | tee -a $out_log_path

echo -e "Running dem_import ...\n" | tee -a $out_log_path
dem_import $in_dem_path $tmp_dem_path $tmp_dem_par_path 0 1 $geoid_dem_path $geoid_dem_par_path 0 | tee -a $out_log_path
echo -e "\nFinish dem_import!" | tee -a $out_log_path

# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
# 重采样 DEM
# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
out_resamp_dem_path=$in_output_folder/$in_produce_name"_resamp.dem"
out_resamp_dem_par_path=$in_output_folder/$in_produce_name"_resamp.dem_par"

echo -e "Running dem_trans ...\n" | tee -a $out_log_path
if [ -f "$out_resamp_dem_par_path" ]; then
    rm -f $out_resamp_dem_par_path
fi
dem_trans $tmp_dem_par_path $tmp_dem_path $out_resamp_dem_par_path $out_resamp_dem_path $in_dem_resamp_factor_y $in_dem_resamp_factor_x - - - - - | tee -a $out_log_path
echo -e "\nFinish dem_trans!" | tee -a $out_log_path

# 可通过 disshd 查看 DEM
# (1) disshd $out_resamp_dem_path 2464 2.778038499e-04 -2.778599818e-04
# (2) 推荐: disdem_par LT1_143448_resamp.dem LT1_143448_resamp.dem_par

# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
# 删除临时文件
# <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
rm -f $tmp_dem_par_path
rm -f $tmp_dem_path
